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ABSTRACT 

Supernova-driven outflows from early galaxies may have had a large impact on the kinetic and chemical 
properties of the intergalactic medium (IGM). We use three-dimensional Monte Carlo cosmological 
realizations of a simple linear peaks model to track the time evolution of such metal-enriched outflows 
and their feedback on galaxy formation. We find that at most 30% of the IGM by volume is enriched to 
values above 10~ 3 Zq in models that only include objects that cool by atomic transitions. The majority 
of enrichment occurs relatively early (5 z 12) and leads to a mass-averaged cosmological metallicity 
between 10 _3 Z Q and 10~ 15 Zq. The inclusion of Population III objects that cool through H2 line 
emission has only a minor impact on these results: increasing the mean metallicity and filling factor 
by at most a factor of 1.4, and moving the dawn of the enrichment epoch to z ~ 14 at the earliest. 
Thus enrichment by outflowing galaxies is likely to have been incomplete and inhomogeneous, biased to 
the areas near the starbursting galaxies themselves. Models with a 10% star formation efficiency can 
satisfactorily reproduce the nearly constant (2 < z < 5, Z rj 3.5 x 10~ 4 Zq) metallicity of the low column 
density Lya forest derived by Songaila (2001), an effect of the decreasing efficiency of metal loss from 
larger galaxies. Finally, we show that IGM enrichment is intimately tied to the ram-pressure stripping 
of baryons from neighboring perturbations. This results in the suppression of at least 20% of the dwarf 
galaxies in the mass range ~ 3 x 10 8 — 3 x 10 9 M & in all models with filling factors greater than 2%, and 
an overall suppression of ~ 50% of dwarf galaxies in the most observationally-favored model. 

Subject headings: intergalactic medium - galaxies: interactions - galaxies: evolution - large-scale 
structure of the universe 



1. INTRODUCTION 

Recent QSO absorption line observations have shown 
that the intergalactic medium (IGM) is polluted with 
heavy elements at intermediate redshifts (Songaila & 
Cowie 1996). From such measurements of column den- 
sity ratios A r civ/A r Hi, Hellsten et al. (1997) and Rauch, 
Hachnelt, & Steinmetz (1997) concluded that typically 
[C/H] ~ —2.5 at z ~ 3, with an order of magnitude dis- 
persion about this mean value. 3 These values, however, 
refer to overdense regions of the universe, traced by Lya 
clouds with column densities in excess of log A/hi = 14.5. 

The presence of metals has more recently been assessed 
in clouds in which log A/hi < 14.0 as reviewed by Pettini 
et al. (2001). At these low optical depths, statistical tech- 
niques to extend the search for highly ionized species such 
as Civ and Ovi must be applied. The results show that 
i) unrecognized weak Civ systems must be present in or- 
der to reproduce the full Civ optical depth (Ellison et al. 
2000), and ii) that metals, as traced by Ovi are present 
in a gas with a density lower than that of the mean IGM 
(Schaye et al. 2000). Very recently, Songaila (2001) has 
been able to trace the IGM metallicity evolution in sys- 
tems with log Civ > 12 and conclude that a minimum 
metallicity Z ss 3.5 x 10~ 4 Z© is already in place at z = 5. 

Although these techniques help to extend QSO absorp- 
tion studies to underdense regions of the IGM, present 
observations are only able to place a lower limit on the to- 
tal volume filling factor of metals. Current measurements, 



combined with numerical simulations, indicate that met- 
als associated with log A/hi ~ 14.2 filaments fill <L 3% of 
intergalactic space, including areas far away from the high 
overdensity peaks where galaxies form (Madau, Ferrara, 
& Rees 2001, hereafter MFR). This suggests that metal 
pollution occurred relatively early, resulting in a more 
uniform distribution and enriching vast regions of inter- 
galactic space. This allows the Lya forest to be hydrody- 
namically 'cold' at low redshifts, as intergalactic baryons 
have time to relax again under the influence of dark-matter 
gravity. Note that the presence of high-redshift metals is 
of great observational importantance, as the measurement 
of metal lines in z ^ 6 quasars may also serve as a probe 
of reionzation (Oh 2002). 

These observations prompted MFR to suggest high- 
redshift (z fts 10) galaxy outflows as a mechanism for IGM 
enrichment. This study could not determine the metal fill- 
ing factor produced in such a scenario, however, as it was 
focused on the evolution of typical objects at a single mass- 
scale and formation redshift. Similar outflow models have 
been proposed by Nath & Chiba (1995) and Scannapieco 
& Broadhurst (2001, hereafter SB), but primarily moti- 
vated by the chemical and thermal properties of the X-ray 
emitting gas in galaxy clusters. While the latter of these 
studies included a range of galaxy masses and was able to 
make some estimates as to the total filling factor of metals, 
these results were fairly crude as the study was focused on 
the properties of individual galaxies. Aguirre et al. (2001a) 
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and Aguirre et al. (2001b) studied IGM metal enrichment 
by superimposing an outflow model on numerical simu- 
lations that did not include SN-driven winds, but were 
only able to constrain the contribution from late-forming 
(z < 6) and relatively large (M > 10 8 5 M Q ) objects. Cen 
& Ostrikcr (1999) studied metal enrichment in even lower 
resolution smoothed particle hydrodynamic (SPH) simu- 
lations with a dark matter particle mass of 8.6 x 10 8 M Q . 
Gnedin & Ostriker (1997) studied the relationship between 
rcionization and early metal enrichment in high-resolution 
simulations, but did not adequately follow supernova ex- 
plosions. Finally, Thacker, Scannapieco, & Davis (2002) 
were able to estimate the filling factor of outflows at z > 4 
purely in the context of high-resolution SPH simulations 
with a dark matter particle mass of 2.5 x 1O 6 M0, but were 
not able to examine its dependence on model parameters 
due to the high computational cost of this approach. 

Early-enrichment scenarios also have important impli- 
cations for the thermal and velocity structure of the IGM, 
as first studied in Tegmark, Silk, & Evrard (1993) and 
Voit (1996) (see also Cen and Bryan 2000). The resulting 
feedback on galaxy formation was first examined in Scan- 
napieco, Ferrara, & Broadhurst (2000), SB, and Scanna- 
pieco, Thacker, & Davis (2001). The nature of this effect 
is twofold: an impinging wind may shock-heat the gas 
of a nearby perturbation to above the virial temperature, 
thereby mechanically evaporating the gas, or the gas may 
be accelerated to above the escape velocity and stripped 
from the perturbation entirely. The latter channel is con- 
siderably more effective, because shock-heated clouds that 
are too large to be stripped are able to radiatively cool 
within a sound crossing time, thus limiting evaporation. 
Note that this type of feedback is fundamentally differ- 
ent from the one commonly adopted in galaxy formation 
models, in which hot gas is produced by supernovae in the 
parent galaxy. 

In this paper we return to the issues of enrichment and 
feedback, adopting a more complete approach that com- 
bines the detailed modeling of a typical object as in MFR, 
with the more general spatially dependent modeling de- 
scribed in SB. In this way we are able place constraints on 
the overall metal filling factor produced as well as inves- 
tigate the link between cosmic metal enrichment and the 
feedback from outflows on galaxy formation. 

The structure of the paper is as follows. In §2 and §3 we 
describe our semi-analytical simulations of galaxy forma- 
tion with feedback and IGM enrichment. In §4 we sum- 
marize the results of these simulations and the constraints 
they place of the fraction of the universe impacted by out- 
flows; conclusions are given in §5. 

2. A LINEAR PEAKS MODEL OF GALAXY FORMATION 

In order to determine the distribution of outflows as a 
function of cosmic time, we use the linear peaks model de- 
scribed in detail in SB. Note that it is important not only 
to have a measure of the overall number density of such ob- 
jects, but also of their spatial distribution, as high-redshift 
galaxies are expected to be highly clustered both from the- 
ory (e.g., Kaiser 1984) and observations (Giavalisco et al. 
1998). 

Using a standard fit to the Cold Dark Matter (CDM) 
power-spectrum (Bardeen et al. 1986), we construct a 256 3 
linear density field spanning a (4 h^ 1 Mpc) 3 cubic comov- 



ing volume, where h is the Hubble constant in units of 100 
km s _1 Mpc . Based mainly on the latest measurements 
of Cosmic Microwave Background (CMB) anisotropies (eg. 
Balbi et al. 2000; Netterfield et al. 2001; Pryke et al. 2002) 
and the abundance of galaxy clusters (Viana & Liddlc 
1996), we focus our attention on a cosmological model with 
parameters h — 0.65, £1m = 0.35, ft a = 0.65, = 0.05, 
(is = 0.87, and n = 1, where Om, ^a, and Sib are the 
total matter, vacuum, and baryonic densities in units of 
the critical density, as is the mass variance of linear fluc- 
tuations on the 8/i -1 Mpc scale, and n is the tilt of the 
primordial power spectrum. 

This linear density field is convolved with spherical 'top- 
hat' window functions corresponding to nine different to- 
tal masses, spaced in equal logarithmic intervals from 
3.0 x 10 7 M to 4.3 x 10 11 M , and spanning the interest- 
ing range from objects that lie close to the lower limit set 
by photoionization and molecular cooling (e.g., Barkana 
& Loeb 1999; Ciardi, Ferrara, & Abel 2000) to the most 
massive galaxies that host outflows. We assume that the 
formation of Population III (Pop III) objects, defined as 
halos with virial temperatures below 10 4 K, is completely 
suppressed by photodissociation of hydrogen molecules by 
UV radiation produced by nearby objects, and we study 
the impact of this assumption in further detail below. 

Using the elliptical collapse model of Sheth, Mo, & 
Tormcn (1999), we identify all 'virialized' peaks in the 
overdensity field, arrange them in order of decreasing col- 
lapse redshift, and exclude all unphysical objects collaps- 
ing within more massive, already virialized halos. After 
collapse, we account for the finite gas cooling time using 
a simple inside-out collapse model (White & Frenk 1991; 
Somcrville 1997). The cooling gas initially relaxes to an 
isothermal distribution at the virial temperature T v j r in 
a Navarro, Frenk, & White (1997, hereafter NFW) dark 
matter halo with concentration parameter c = 5, and with 
a uniform metallicity Z calculated as described in the next 
Section. In this model, the gas within a radius r coo i cools 
due to radiative losses that account for metallicity as tab- 
ulated by Sutherland and Dopita (1993), and the gas out- 
side this radius stays at the virial temperature of the halo, 
with r coo i moving outwards with time. This model is de- 
scribed in further detail in SB. 

When the total mass contained within r coo i equals the 
object's baryonic mass, a new galaxy is assumed to form 
with a gas mass Mb — (Slb/SlM)M. In the smaller ha- 
los at high redshift having logT v i r < 5.7 (i.e. masses 
M < 2xl0 10 [(l+z)/10]" 3 / 2 M ), rapid cooling by atomic 
hydrogen and helium occurs on timescales much shorter 
than the gas free-fall time, and infalling gas collapses to 
the center at the free-fall rate rather than coming to hy- 
drostatic equilibrium (MFR). The supply of cold gas for 
star formation is then only limited by the infall rate. 

We assume a Salpeter initial mass function (IMF) with 
upper and lower mass cut-offs equal to M u — 120 M Q and 
Mi =0.1 M©, respectively. In this case, one SN occurs for 
every v^ 1 k, 136 M Q of stars formed, releasing an energy 
of Eo = 10 51 erg. This mechanical energy is injected by 
SNe after a few times 10 7 yr: at this stage SN-driven bub- 
bles propagate into the halo quenching further star forma- 
tion, and the conversion of cold gas into stars is limited 
by the increasing fractional volume occupied by SN rem- 
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nants. To be conservative, we do not consider the possibil- 
ity that very massive (« 300M Q ), metal-free stars might 
contribute both to mctallicity and energy input. Their 
effects have been discussed in detail by Schneider et al. 
(2002) and Schneider, Guetta & Ferrara (2002) to which 
we refer to the reader for an extensive description. 

Before SN feedback occurs some fraction /* of the gas 
will be able to cool, fragment, and form stars. As our for- 
malism does not include local feedback effects, this star 
formation efficiency must be considered as a free parame- 
ter of the model. Finally, a fraction /„, of this energy will 
be channeled at a constant rate into a galaxy outflow over 
a timescale span of ioB = 33 Myr, ejecting gas into the 
IGM. 

As we will show below (Fig. 1), the most efficient IGM 
pollutors are objects with masses of a few times 10 8 M Q , 
for which the fraction of gas that can cool in a free-fall 
time is essentially unity (MFR). Hence, this gas is readily 
available to be transformed into stars on short timescales. 
This justifies the prompt star formation (starburst mode) 
approximation we have adopted, which therefore should 
be appropriate to the aims of this study. 

3. MODELING GALAXY OUTFLOWS 

The outflows arc modeled as spherical shells using a 
method that is based on the approach described in SB, 
but with several important refinements taken from MFR, 
Ferrara, Pettini & Shchekinov (2000), and Mori, Ferrara 
& Madau (2001; hereafter MFM). An outflow is driven 
out of the galaxy by internal pressure and decelerated by 
inertia and the gravitational pull of the dark matter halo, 
both estimated in the thin shell approximation (Ostriker 

6 McKee 1988; Tegmark, Silk, & Evrard 1993). 

The expansion of the shell, whose radius is denoted by 
R s , is driven by the internal energy, Et, of the hot bub- 
ble gas. The pressure of such a gas (with adiabatic index 

7 = 5/3) is therefore = Ei,/2ttR^. Momentum and 
energy conservation yield the relevant evolutionary equa- 
tions: 



R s — ' - b „ — tt(Rs ~ HR S ) 2 - Q. M — — - 

pR s R s 2 

E b = L{t) - 4irR 2 s R s P h - L c , 



(1) 



where the dots represent time derivatives, the subscripts 
s and b indicate shell and bubble quantities respectively, 
g s = GM(R S )/R 2 , and p is the density of the ambient 
medium, taken to be the halo gas density within the virial 
radius and the mean IGM background density outside the 
virial radius. These equations reduce to those given in 
MFR in the regime in which the Hubble expansion is neg- 
ligible and reduce to those given in SB if the NFW profile 
is replaced by a point mass and external pressure is ne- 
glected. 

The cooling rate, L c , is assumed here to be dominated 
by inverse Compton cooling off CMB photons (Ikeuchi & 
Ostriker 1986), as gas radiative processes are much less effi- 
cient in the low density 10 5 K < T < 10 8 K gas that drives 
the outflows. This approximation is especially appropri- 
ate as the combined cooling processes produce variations 
of less than a few percent on the final size of the bubble 
(see Fig. 6 of MFR). 



The mechanical luminosity of SNe is given by 
L(t) = (f w E )^Q(t OB - t) 

tOB 

= 8.26 x 10 33 9(t O B - t) ( )f w f*M ergs" 1 . (2) 

This assumption of a constant luminosity over the burst is 
most accurate for the larger galaxies in our simulations, in 
which the stochastic variations of L(t) become smaller due 
to the larger number of SNe. We constrain f w by combin- 
ing the overall efficiency of 30% derived for the 2 x 10 8 M Q 
object simulated by MFM with the mass scaling derived 
in Ferrara, Pettini & Shchekinov (2000), which was ob- 
tained by determining the fraction of starburst sites that 
can produce a blowout in a galaxy of a given mass. Thus, 
we choose f w (M) = 0.35 B {M)/8 B (M = 2xl0 8 M Q ) where 



5 B (M) 
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[1.0- 0.165 ln(100)] 100 N t ~ 



N t < 1 
1 < N t < 100 

100 < N t 

(3) 

where N t = 1.7 x 10 7 {n b /D, M )M /M Q is a dimensionless 
parameter that scales according to the overall number of 
SNe produced in a starburst, divided by the efficiency /*. 

Within the virial radius a fixed fraction f m = 0.5 of the 
gas is swept into the shell, a value taken from the numer- 
ical simulations described in MFM. In those experiments 
it is seen that after blow-out, half of the initial mass con- 
tained in the virial radius re-collapsed to the center as a 
result of the multiple shell-shell interactions leading to the 
formation of cold sheets. In this case the halo gas is as- 
sumed to virialize to an isothermal distribution and settle 
down to a density profile 



lnp(r) = ln(/ m( o ) 



m-vl(r)}, (4) 



(Makino, Sasaki, & Suto 1998), where the central, pre- 
burst gas density po is determined by the condition that 
the total baryonic mass fraction within the virial radius is 
equal to the cosmic average, yielding p$ = 11052 / o cr i t 5!(,- 

Outside the virial radius the shells expand into the Hub- 
ble flow, sweeping up all of the baryons in their path. Fi- 
nally, when outflows slow down to the point that they are 
no longer supersonic, our approximations break down, and 
the shell is possibly fragmented by random motions. At 
this point we let the bubble expand with the Hubble flow. 

To calculate the cooling time of forming galaxies (see 
§2), we use a simple estimate of the metallicity of col- 
lapsed halos. The metals in each outflow are assumed to 
be evenly distributed within its radius, adopting an aver- 
age yield from each supernova of 2 M Q (e.g. Nagataki & 
Sato 1997), half of which is deposited in the outflow and 
half of which remains in the galaxy. Each collapsing ob- 
ject is then assigned a mass in metals Mz, taken to be zero 
initially and modified by each outflow passing within its 
collapse radius, r co n. We assume that the fraction of met- 
als which fall into the collapsing halo from a given outflow 
is equal to the volume fraction of the outflowing bubble 
which falls within the collapsing sphere. In this case, for 
each overlapping bubble, Mz is updated to 
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where K> V criap is the volume of intersection between the 
outflow and the collapsing sphere. By dividing this mass 
by the total baryonic mass of the galaxy we can compute 
the initial metallicity of the object, and thus the delay be- 
tween collapse and star formation due to cooling. Note 
however, that this cooling time is only significant at the 
highest mass scales, ^ 1O 1O M0. 

Our model is also able to account for the inhibition of 
surrounding low-mass galaxy formation (i.e. negative feed- 
back). Scannapieco, Ferrara, & Broadhurst (2000) showed 
that the most important such mechanism is 'baryonic 
stripping,' whereby high-redshift galaxy outflows strip the 
gas out of nearby overdense regions that would have other- 
wise become low-mass galaxies. Whenever a shock moves 
through the center an overdense region that has not yet 
virialized, we apply a simple check to determine if such 
stripping has occurred. 

In each case we estimate the comoving radius of the 
overdense region as r/w (1 + ^nl) -1 / 3 where i5nl is the non- 
linear overdensity of the region, which estimated from the 
spherical collapse model as 



1 + 6nl = 



2 (1 - cos* 



(6) 



where the collapse parameter 8 is given by (9 — 
sin0) 2 / 3 7r~ 2 / 3 = D(z cmss )/D(z c ) where z cross is the red- 
shift at which the shock moves through the center of the 
objects, z c is its collapse redshift, computed from the 
model of Sheth & Tormen (2001), and D(z) is the linear 
growth factor as a function of redshift. We then inhibit 
galaxy formation if the shock has sufficient momentum to 
accelerate the gas to its escape velocity: 



> M p V e 



(7) 



where u> is the solid angle of the shell subtended by the 
perturbation, M s is the baryonic mass swept up by the 
shock, including both the IGM and gas ejected from the 
host objects, v s the velocity of the shock, and M p and 
v e = 2GM P (1 + ^nl) 1 / 3 / r M p are the gas mass and escape 
velocity of the perturbation. 

In SB we simply estimated the filling factor of outflows 
by summing over the volumes contained within all galaxy 
outflows, regardless of their mass scales or spatial distri- 
bution. As this investigation is focused on determining 
the overall state of the IGM, here we adopt a more accu- 
rate approach in which we construct a 256 3 grid, search 
through the expanding bubbles, and count all grid points 
within them as a function of redshift . 

Our approach also differs from that in SB in our treat- 
ment of outflows expanding from host galaxies that merge. 
While our previous estimates only counted outflows from 
unmcrged host galaxies, in this work we instead take the 
wind to remain fixed at the comoving radius found at the 
time of the merger. This is because, while it is clear that 
our approximation of a spherical pressure-driven wind is 
broken when a merging galaxy passes through this shell, it 
is unlikely that such a merger would be able to absorb all 
the gas ejected, somehow erasing the wind from existence. 

It is important to point out that although our estimate 
of the filling factor takes into account the spatial orien- 
tation of the bubbles, our formalism fails to capture the 



interaction between expanding shock fronts. In the case in 
which two outflows are expanding in opposite directions, 
our approach allows them to pass through each other un- 
hindered, leading to a slight overestimate of the volume 
impacted. Similarly, in the case in which an outflow ex- 
pands into a bubble evacuated by a previous wind, it is 
nevertheless modeled as if it were expanding into the Hub- 
ble flow, leading to a slight underestimate. 



4. RESULTS 



4.1. 



Cosmic Metal Filling Factor 

In the left panels of Figure 01 we show the derived fill- 
ing factor of outflows. Here, the bend at z ~ 13 is due 
to our finite mass resolution. At this redshift, the small- 
est mass scale drops below the virial temperature limit of 
10 4 K, causing the effective resolution of the simulation to 
jump to the next mass scale. In reality objects between 
these two mass scales are formed, and the metals from 
these objects would smooth the evolution of the filling 
factor. Finally, because our method becomes inaccurate 
when perturbations are overly nonlinear, we restrict this 
figure to redshifts above 3. 
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Fig. 1. — Left panels: Cosmic metal filling factor from galaxy 
outflows. In each panel solid (dashed) lines show the filling factor 
in runs with (without) the inclusion of baryonic stripping feedback 
according to eq. (pj) . Each panel is labeled by the overall star forma- 
tion coefficient and in the /* =0.1 case the dotted lines represent 
models with no feedback, but outflows are only allowed from objects 
at or above a fixed mass limit. From top to bottom these limits are 
1.0 x 1O 8 M , 3.3 x 1O 8 M , 1.1 x 1O 9 M , and 3.6 x 1O 9 M . Right 
panels: Mass-averaged IGM metallicity as a function of redshift. 
Curves are as in the left panels. 

The parameters in these models are /* = 0.5,0.1,0.01, 
L,(M) = 0.38 B (M)/S B (2 x 10 8 M q ) with 6 B as in eq. 
(§), v- 1 = 136M , t B = 33 Myr, and f m = 0.5. The 
first three of these parameters appear in combination both 
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in the expression for the overall strength of the winds 
(oc f-kfwv), and the metallicity (oc f+v). On the other 
hand, ioB has almost no effect on our results as the rele- 
vant times for star formation are small compared to struc- 
ture formation times scales. Thus we can provide a con- 
servative estimate of the model uncertainties introduced 
by these parameters by simply considering a wide range of 
star formation efficiencies, and applying a linear shift in 
the final metallicity to estimate the effect of varying f w . 
Finally, while the mass loading parameter f m has little ef- 
fect on the overall filling factor, it is important for galaxy 
feedback, and we consider its impact in detail §4.3. 

The most obvious, yet perhaps most important feature 
of Fig. 1 is that the filling factor is always substantially 
less than unity, ranging from 3% to 30% at z — 3. Note 
that these values are consistent with the 20% enrichment 
at z = 4 found in numerical simulations by Thacker, Scan- 
napieco, & Davis (2002), using a model similar to our 
/* = 0.1 case. The fact that IGM enrichment is inho- 
mogeneous even in the maximal case in which 50% of all 
baryons in collapsed objects are taken to form stars, how- 
ever, leads us to an important conclusion: starburst driven 
outflows, while an effective source of metals in overdense 
regions (SB), are not able to enrich the IGM in its en- 
tirety. This is true even in the ACDM model considered 
in our simulations, in which dwarf galaxies are formed at 
very high redshifts, and the baryonic/dark matter ratio is 
relatively high, resulting in a large number of stars. 

The details of our results depend sensitively on the mini- 
mum mass scale of the galaxies in our simulation, however, 
which is set by our minimum virial temperature of 10 4 K. 
In the central panel of this figure, we plot a series of mod- 
els in which no feedback as per eq. (|?|) is imposed, but 
instead we allow outflows only from objects above a fixed 
mass scale. Both the redshift at which outflows begin to 
become important and their overall filling factor depends 
closely on this mass. Thus, while in the run with /* = 0.1, 
outflowing bubbles fill 1% of the volume at redshift ;$ 12 
and reach a final filling factor of 16%, excluding all ob- 
jects with masses below 1.1 x 10 9 M Q shifts these values to 
z « 8 and 6% respectively. Note that this lower resolution 
is similar to that adopted by Aguirre et al. (2001a) and 
approximately equal to the mass of a single dark matter 
particle in the simulations by Cen and Ostriker (1999). 

In spite of the sensitivity of metal enrichment to low- 
mass objects, its overall dependence on baryonic stripping 
feedback is weak, as can be seen by comparing the solid 
lines in which equation (7j) has been imposed with the 
dashed lines in which such feedback from outflows is ne- 
glected. The shape and final value of the filling factor are 
extremely similar between such models for all values of 
/*, becoming indistinguishable in many cases. This is be- 
cause baryonic stripping can only occur in a perturbation 
that is sufficiently nearby and late-collapsing. Then the 
shock velocity, R s , is large and the overdense region oc- 
cupies a large solid angle, w, when the outflow reaches it. 
Thus the perturbations succumbing to baryonic stripping 
correspond to late-forming galaxies in the most heavily 
populated regions of space, which have little effect on the 
overall filling factor. 

The higher bias of suppressed objects can also be seen 
by comparing the evolution of the filling factor with the 



overall mass-averaged IGM metallicity, plotted in the right 
panels of Figure [y. In these panels, the differences be- 
tween the models with and without suppression are much 
more pronounced. The difference is most apparent in the 
/* = 0.5 case, in which the wind velocities are the high- 
est, and thus the suppression of neighbors is most severe. 
In this case at z = 3 the overall metallicities differ by a 
factor of 1.5 while the difference in volume filling factor is 
less than a factor of 1.15. 

Note that the mass-averaged metallicity scales almost 
linearly with as this parameter controls the number 
of stars formed in each galaxy, and hence the number of 
supernovae and mass of ejected metals. We find that at 
z = 3, Z m 0.1/*, where this relation depends on the as- 
sumed yield (2M Q per SN, 1/2 ejected), the gas ejected 
fraction (50%), and the minimum mass scale in the sim- 
ulation. This mass dependence, while sensitive, is more 
limited than that of the overall filling factor, as can be 
seen by comparing the /* = 0.1 model with the series of 
models with a threshold mass imposed, plotted in the cen- 
ter right panel. Thus the model with T v i r > 10 4 K (solid 
line) and that with a fixed threshold of 1.1 x 10 9 M Q have 
z = 3 metallicities of 0.010 Zq and 0.007 Z Q respectively. 
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Fig. 2. — Cosmic metal filling factor from galaxy outflows, as a 
function of threshold metallicity. In each panel the solid lines show 
the total filling factor of metals, which is indistinguishable from the 
filling factor of gas with Z > 10~ 3 Zq. The dotted lines, from top 
to bottom, show the filling factor of gas enriched to values above 
(10" 2 , 1CT 1 , 1)Zq, respectively. 

In Figure ^ we examine the distribution of metals in 
more detail, by computing the filling factor of gas above 
a fixed threshold metallicity. As the metallicities of indi- 
vidual outflows always exceed 10~ 3 Zq in all models, the 
filling factor at this threshold metallicity is indistinguish- 
able from the total filling factor as shown in Figure [y. 

Choosing a higher threshold value, however, leads to 
more model-dependent results. From this point of view, 
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the /* = 0.5 and /* = 0.1 models are quite similar, and 
in both cases approximately 50% of the enriched gas has 
a metallicity greater than 1Q~ 2 Zq, 10% of the enriched 
gas exceeds 10 _1 Z Q , and 1% of the enriched IGM exceeds 
solar metallicity. On the other hand, the distribution of 
metals is much more clumpy in the weak-star formation 
(/* = 0.01) case. Here, while about 50% of the enriched 
IGM exceeds 10~ 2 Z Q , the volume fraction of enriched gas 
with metallicities above 10 Zq and Zq are quite similar, 
with such high-metallicity gas occupying approximately 
6% and 3% of the enriched volume respectively. 

Our estimates of the IGM metallicity invite comparison 
with that recently observed in the low-column density Lya 
forest. Songaila (2001) derived the redshift evolution of 
ficrv an d figuv f° r a large number of systems (367 and 109 
for Civ and Siiv ) with H i column densities 10 12 cm~ 2 < 
JVhi < 10 15 cm -2 . The above quantities can then be trans- 
lated, modulo the assumption of the ionization fraction of 
such species taken to be equal to 0.5, into a metallicity 
estimate. The resulting metallicity evolution is shown by 
the points in Fig. ^, with the error bars corresponding to a 
90% confidence limit. Songaila points out two striking fea- 
tures of IGM enrichment: (i) the existence of a minimum 
metallicity Z sa 10~ 3A5 Z & , already in place at z = 5, and 
(ii) a constancy of metal abundances through the studied 
redshift range z = 2 — > 5. A direct comparison of our 
simulations with these results is possible provided that we 
filter our metallicity estimates with the overdensities sam- 
pled by Songaila's experiment. Ricotti, Gnedin & Shull 
(2000) showed that a clear correlation between H i column 
and gas overdensity exists, with Pbl~p~b — 0.8 iVjjj 7 13 ; hence, 
the 10 12 cm' 2 < Nm < 10 15 cm~ 2 range can be translated 
into the overdensity range 0.16 ,$ Pb/p~b ~ 20. 
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Fig. 3. — Metallicity evolution for low density (10 12 cm 2 < 
N m < 10 15 cm~ 2 ) IGM. Points are the data of Songaila (2001); 
dotted, solid, and dashed lines are the predicted values for = 
0.5,0.1,0.01 respectively. 



Adopting the smallest filter scale and associating lin- 
ear and nonlinear overdensities according to eq. (||), 
we have searched through our simulations and com- 
puted a mass averaged metallicity of cells that are not 
surrounded by a collapsed perturbation and in which 
0.8 Pbl~Pb ~ 20, roughly corresponding QSO absorption 
systems with 10 13 cm -2 < A ffl < 10 15 cm~ 2 . We adopt 
this lower bound as our mapping between linear and non- 
linear overdensities becomes unreliable for smaller values 
of 5. These results are superposed on the data points in 
Fig. H for the usual three values of Although no fit- 
ting parameter has been adjusted to produce these curves, 
the resemblance to the data, both in terms of amplitude 
and flat shape, is noticeable. It is only at redshifts above 
z — 8 that the curve starts to bend downwards. The phys- 
ical interpretation for the existence of this metallicity floor 
is simple: metal enrichment is dominated by low-mass 
(10 s — 10 9 M©) objects for which metal ejection/transport 
occurs with highest efficiency (see discussion in Ferrara, 
Pettini & Shchekinov and eq. 0). As the nonlinear mass 
scale increases, larger but intrinsically less efficient galax- 
ies become the most numerous population. However, they 
contribute only marginally to metal enrichment of the Lya 
forest as synthesized metals are predominantly trapped by 
their potential wells. This explains both the shape and 
amplitude of the data, provided /* « 0.1, a value that is 
broadly consistent with other arguments based on cosmic 
star formation rates (Ciardi et al. 2000; Thacker, Scanna- 
pieco, & Davis 2002; Barkana 2002). 



4.2. Population III Objects 

In order to study the impact of Pop III objects on our 
results, we have conducted a number of comparison simu- 
lations in which we allow objects with virial temperatures 
< 10 4 K to form, but with a lower overall star formation 
efficiency. In these cases we took f+popin = 0.1/* in ob- 
jects with T < 10 4 K. This is physically motivated by the 
lower efficiency of molecular hydrogen cooling, which re- 
duces the number of cooled baryons available to form stars. 
For each value of /* we conducted three such comparison 
runs: one with the same overall resolution as the fiducial 
case (a 4 h~ Y Mpc comoving box with objects with nine 
mass scales ranging from 3.0 x 10 7 M Q to 4.3 x 10 11 M Q ), 
one in a 2.75 h~ x Mpc comoving box with the nine mass 
scales now ranging from 1.0 x 10 8 M Q to 1.3 x 10 11 M , and 
one in a 1.75 Mpc comoving box with masses rang- 
ing from 3.0 x 10 6 M Q to 4.3 x 10 11 M . Note that in the 
last of these simulations the smallest objects have virial 
temperatures < 5, 000 K at redshifts below 20, and are 
thus not only well below the 10 4 K atomic cooling limit, 
but also below the 5, 000 K limit for "efficient" molecular 
cooling at which the cooling time exceeds the free-fall time 
(MFR). 

The filling factor and overall metallicities in the com- 
parison runs are shown in Figure ^. Here we restrict our 
analysis to redshifts above 5, so that our linear approach 
remains accurate at the smallest mass scales. In this figure 
we see that even in the case with the highest resolution, 
Pop III objects are only able to increase the metal filling 
factor and metallicity of the IGM to values that are a fac- 
tor of 1.4 times greater than those found in the standard 
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runs. Similarly, these objects have only a small impact on 
the redshift at which outflows fill 1% of the volume, push- 
ing this value from 12 to 14 in the /* = 0.1 case. From 
this comparison, it is clear that while our results are sensi- 
tive to imposing a mass limit above T V ir = 10 4 K, they are 
relatively insensitive to the inclusion of Pop III objects, 
which rely on H2 line cooling for collapse. 
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Fig. 4. — Impact of Pop III objects. In each panel the solid lines 
show the filling factor and total metallicity in runs with a minimum 
mass scale of 3.0 X 10 7 Mq and a minimum T v j r = 10 4 K as in Fig. 
1. The dotted lines are taken from runs which include star forma- 
tion in objects T v i r < 10 4 K as described in the text. From top to 
bottom the minimum mass scale in these simulations is 3.0 X 10 6 
Mq, 1.0 X 10 7 Mq, and 3.0 X 10 7 Mq, corresponding to virial tem- 
peratures of 5,000, 12,000, and 26,000 K at 2 = 20 and 2,800, 6,000, 
and 13,000 K at z = 10. 



4.3. Feedback 

In Figure || we explore in greater detail the relationship 
between the cosmic filling factor and baryonic stripping. 
While the suppression of galaxy formation has little im- 
pact on the cosmic metal filling factor, these quantities 
are clearly related, as both are dependent on the number 
density and strength of galaxy outflows. In this plot we 
quantify baryonic stripping by counting all galaxies with 
masses above 3.2 x 10 8 M Q , excluding the objects at the 
smallest two mass scales, as at low redshifts their virial 
temperatures fall below the 10 4 K limit. In this case, the 
majority of suppressed objects are slightly larger dwarf 
galaxies with masses of a few times 10 9 M Q (SB; Scanna- 
pieco, Thacker, & Davis 2001). 

In the upper panel of this figure, we see that while the 
details of outflow generation and propagation introduce 
a large scatter, widespread IGM enrichment is accompa- 
nied by a significant level of baryonic stripping in all cases. 
Thus even given the wide range of /+ values considered, all 
models and redshifts in which 2% of the IGM is enriched 



show a suppression of galaxies £ 20%. In the /* = 0.1 
case, which is most consistent with QSO observations, ap- 
proximately half of the objects with masses at or above 
3.2 x 1O 8 M0 are suppressed by this mechanism by z = 3. 

These results are sensitive however, to the mass-loading 
parameter f m , as the suppression of galaxy formation is 
primarily due to the momentum carried by the winds (see 
eg. Scannapieco, Ferrara, & Broadhurst 2000; Scanna- 
pieco, Thacker, & Davis 2001), which scales as E^ Ml . 
Thus increasing the fraction of gas swept into the shells in- 
creases galaxy suppression in models with the same kinetic 
energy input. In the lower panel of Figure [5], we compare 
our fiducial /* = 0.1, f m = 0.5 model with enhanced (de- 
pressed) feedback models in which f m = 1.0 (/,„ = 0.25). 
Although changes in /+ result in a simultaneous increase 
of both the filling factor and suppression factor of dwarf 
galaxies, changing f m shifts the suppression factor while 
having little effect on the overall filling factor. Thus the 
filling factors in these models are all between 14% to 18%, 
while the fraction of suppressed dwarf galaxies varies from 
35% to 60%. 
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Fig. 5. — Top: Fraction of galaxies suppressed by baryonic strip- 
ping as a function of the overall metal filling factor. The dotted, 
solid, and dashed curves correspond to models in which /* is set 
to 0.5, 0.1, and 0.01 respectively. Bottom: Impact of f m on the 
suppression of galaxy formation. From top to bottom f m = 1.0, 0.5, 
and 0.25, with /* = 0.1 in all cases. 

As a general remark about this type of feedback, we 
note that the gas is swept from the potential well is un- 
likely to be reaccreted by this object. Instead, this gas be 
available to other collapsing objects, which therefore may 
be characterized by higher baryon-to-dark matter ratios. 
Clearly, in order to draw more quantitative conclusions a 
dedicated hydrodynamical study of the details of baryonic 
stripping is required. 
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5. CONCLUSIONS 

In this work we have studied the metal enrichment of the 
IGM by outflows in a ACDM model of structure formation 
and its feedback on the formation of galaxies. Adopting 
a linear peaks model of the spatial distribution of form- 
ing objects, and a detailed one-dimensional model of wind 
propagation, we have determined the overall filling fac- 
tor as a function of redshift and its relationship with the 
baryonic stripping of protogalaxies. 

While the star formation efficiency of high-redshift 
galaxies is largely unknown, we are nevertheless able to 
place useful constraints on the filling factor, enrichment 
redshift, and overall mass- averaged metallicity in such 
models. Choosing a range of star formation efficien- 
cies between /* = 0.50 and /* = 0.01, we find that at 
least 3% and at most 30% of the IGM is enriched to a 
level exceeding 10~ 3 Zq by redshift z = 3. In all cases, 
the majority of this enrichment occurs relatively early, 
5 ;$ z ^ 12, and leads to mass-averaged cosmic metallic- 
ities that range from 0.001Z Q to O.O5Z , for star forma- 
tion efficiencies 0.01 < /* < 0.5, respectively. The mass- 
averaged metallicity scales roughly linearly with this quan- 
tity: Z«O.1/*Z , 

Our model can satisfactorily reproduce the constant 
(Z « 3.5 x 10 _4 Z Q ) metal enrichment of the low column 
density Lya forest up to z = 5 derived by Songaila (2001), 
which is likely to be caused by the decreasing efficiency of 
metal loss from larger galaxies. This comparison strongly 
favors star formation efficiencies in a narrow range around 



10%, essentially excluding the /* = 0.5 and /* = 0.01 
models. As the formation of stars in Pop III objects is rel- 
atively inefficient, the inclusion of these objects has only 
a secondary effect on our results: increasing the mass- 
averaged metallicity and filling factors by at most a factor 
of 1.4, and moving the dawn of the enrichment epoch to 
z w 14 at the earliest. 

While all the models studied display suppression of 
galaxy formation due to outflows ram-pressure stripping 
the gas out of pre-virialized protogalaxies, this mechanism 
has only a minor impact on the overall filling factor as it 
occurs only in the densest and most polluted regions of 
space. Nevertheless, after fixing f m , a general relation- 
ship between the filling factor and the suppression factor 
of galaxies exists, at all /+ values. All models and redshifts 
at which 2% of the IGM is enriched show a greater than 
20% suppression of galaxies. In the case that is most con- 
sistent with QSO observations, in fact, half the galaxies 
are suppressed due to baryonic stripping. Thus the rela- 
tive quiescence of the Lya forest at lower redshifts is likely 
to belie a violent epoch of early outflows and enrichment. 
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